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Abstract 

The phase-integral approximation devised by Froman and Froman, is used for computing cosmological perturbations in the 
power-law inflationary model. The phase-integral formulas for the scalar and tensor power spectra are explicitly obtained up 
to ninth-order of the phase-integral approximation. We show that, the phase-integral approximation exactly reproduces the 
shape of the power spectra for scalar and tensor perturbations as well as the spectral indices. We compare the accuracy of the 
phase-integral approximation with the results for the power spectrum obtained with the slow-roll and uniform approximation 
j - ■ methods. 
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I. INTRODUCTION 



The results reported by WMAP favor inflation l] over other cosmological scenarios. The data is consistent with 
a flat universe and with an almost scale invariant spectrum for the primordial perturbations. The spectrum of the 
perturbations generated during inflation depends on the model, therefore it is important to predict the power spectrum 
of the cosmological perturbations for a variety of inflationary scenarios. In general, most of the inflationary models 
are not analytically solvable and approximate or numerical methods are mandatory. Traditionally, the method of 
approximation applied in inflationary cosmology is the slow-roll approximation 2]. Recently, some authors have 
applied semiclassical methods, such as the WKB method with the Langer modification [1, 13, [f|, and the method 
of uniform approximation [(J 0, Q. ^ n the present article we propose an alternative method of approximation for 
the study of cosmological perturbations during inflation, this method is based on the phase integral approximation 
0, El , which has been successfully a ppl ied in different problems in quantum mechanics [§| , and in the study of 
quasinormal modes in black hole physics jl2l . [f3j . 

The Friedmann-Robertson- Walker line element for a spatially flat universe can be written as 



ds 2 = -dt 2 



\t) [dr 2 



\de 2 



sin 



edef)} , 



(i) 



where a is the scale factor. In order to study cosmological perturbations we consider perturbations of the spatially-flat 
Friedmann-Robertson- Walker universe (Q]). 

Using a gauge invariant treatment of linearized fluctuations in the metric and field equations [14| ,we have that the 
power spectra involves the computation of the two-point function 



\u(r),x),u(r),x + r)\0) = 



dk sinkr 
k kr 



Pu(v,k). 



(2) 



The scalar density perturbations are described by the function Uk — o-^/ij) 1 , where $ is a gauge-invariant variable 
corresponding to the Newtonian potential, and 4> is the scalar field. The equations of motion for the perturbation Uk 
in a universe dominated by a scalar field <f> are given by 



(3) 



where zs — a<fi' /7i, 7i — a' /a, and the prime indicates derivative with respect to the conformal time rj = J dt/a. 

For the tensor perturbations (gravitational waves) we introduce a function Vk — ah, where h is the amplitude of 
the gravitational wave. The tensor perturbations obey a second order differential equation analogous to Eq. ([3]) 



a 
a 



Vk 



(4) 



Considering the limits k 2 ^S> \z'g/zs\ (short wavelength) and k 2 <C \z'g/zs\ (long wavelength), we have that the 
solutions to Eq. ^ exhibit the following asymptotic behavior 



Uk 



-ikr} 



2k 



(k 2 > \z%/zs\,-kri^KX>) , 



(5) 



A k z s (fc 2 « \4/z s \,-kri^0), 



(6) 



the same asymptotic boundary conditions also hold for tensor perturbations 

Once the mode equations for scalar and tensor perturbations are solved for different momenta k, the power spectra 
for scalar and tensor modes are given by the expression 



Ps{k) = 
P T (k) = 



k 3 

lim — — 

-kr/^O 2-K 1 

k z 

lim — - 



Uk(v) 



a{m) 



(7) 
(8) 
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The spectral indices are denned as [IE 



d\nP s (k) 



d\nP T (k) , . 



Running of the spectral indices is given by the second logarithmic derivative of the power spectra 



dlnns(k) 

ar <*> = ^/H" (12) 



The power spectrum is usually fitted using the ansatz 



k 



nf H +(l/2)a fit ln(fc/fc„) 



P(k) = A /it I — ) , (13) 

where the parameters A fa, Ufa, and a/j* are fitted with the observational data. The spectral index is evaluated at 
the value fc* UM 



ns(K) = 1 + rifu. (14) 

The running is parametrized by a fa. It is not possible to have a constant n and nonzero a given the definitions for 
P(k) and n<?. This inconsistency gives as a result a uncontrolled growth of errors away from the parameter fc». It is 
purpose of the present paper to show how to compute approximate solutions for the scalar and tensor power spectra 
with the help of the phase-integral approximation method. The article is structured as follows: In section II we give an 
introductory review of the phase-integral method and the connection formulas. In Sec. Ill we apply the phase-integral 
approximation to the power-law inflationary model. In Sec. IV we compare the results for the power-spectra obtained 
using the phase-integral approach with those computed with the slow-roll and uniform-approximation methods. 



II. THE PHASE-INTEGRAL METHOD 

Let us consider the differential equation 

d 2 u k 



dz 



+ R(z)u k = 0, (15) 



where R(z) is an analytic function of z. In order to obtain an approximate solution to Eq. (|15[) , we are going to 
use the phase integral method developed by Froman [l7| . The phase integral approximation, generated using a 
non specified base solution Q(z), is a linear combination of the phase integral functions [ToL [T^ |. which exhibit the 
following form 

u fc = ?- 1 / 2 (2)exp[±iw(«)], (16) 

where 



lu(z) = J q(z)dz. (17) 
Substituting (|16[) into (| L 5|) we obtain that the exact phase integrand q(z) must be a solution of the differential equation 



3 



-3/2 



(18) 

For any solution of Eq. (fl8|) the functions ([T6]) , are linearly independent, the linear combination of the functions u k 
represents a local solution. In order to solve the global problem we choose a linear combination of phase integral 
solutions representing the same solution in different regions of the complex plane. This is known as the Stokes 
phenomenon Q. 

If we have a function Q(z) which is an approximate solution of Eq. Q18p . the quantity e , obtained after substituting 
Q(z) into Eq. JT? 



(19) 



is small compared to unity. We take into account the relative small size of eo by considering it proportional to A 2 , 
where A is a small parameter. The parameter e is small when Q(z) is proportional to 1/A and R(z) — Q 2 (z) is 
independent of A, i.e. if R{z) is replaced by Q 2 (z)/X 2 + [R(z) — Q 2 (z)\ in Eq. (jTSJ) . Therefore, instead of considering 
Eq. (|15[) . we deal with the auxiliary differential equation 



d 2 u k 



Q 2 (z) 
A 2 



+ [R.(z)-Q 2 (z)]^u k = 0, 



(20) 



which reduces to Eq. (jl5|) when A = 1. 

Inserting the solutions (|16p into the auxiliary differential equation (|20[) . we obtain the following equation for q{z) 



which is called the auxiliary q equation. After introducing the new variable £, 

f= J Q(z)dz, 

we obtain 



1 - 



q\ 



Q(z) 



e X 2 



_qX_ 

Q(z) 



1/2 l_ 



qX 



1-1/2 



Q(z) 



x 2 = o, 



where eo is defined by Eq. ([19]) . A formal solution of Eq. (|23|) is obtained after the identification 

qX 



E y 2„A 2 



n=0 



Substituting Eq. ([24]) into Eq. ([23]). we obtain 



(21) 



(22) 



(23) 



(24) 



1/2 



E^ 2 " +eoA 2 + E^ 2n j ^[^ Y2nX2n 



-1/2 



(25) 



Using computer manipulation algebra it is straightforward to obtain the coefficients Yi n . The first values are [13, Gil 



Y 
Yi 

n 

la 



32 



' (4 + £ 2) , 

(2e 2 + 6e e 2 + 5e 2 + e 4 ) 



128 



(5e^ + 30e 2 e 2 + 50e o e 2 + 10e e 4 + 28eie 3 + 19e^ + e 6 ) 



(26) 
(27) 

(28) 

(29) 

(30) 
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where e„ is defined as 



1 de "- 1 ^^ /n\ 



Truncating the series (|24|) at n = N with A = 1 we obtain 



N 

q(z) = Q(z)Y,Y2n, (32) 

Substituting (J32} in we have that 



A' 

u{z) = ^2 u 2n {z), (33) 

n=0 

where 



wa»(*) = / r 2n Q(z)^. (34) 

From (|32|) , ([33|l and (fT6|) we obtain a phase integral approximation of order 2N + 1 generated with the help of the 
base function Q(z). 

The base function Q(z) is not specified and its selection depends on the problem in question. In many cases, it is 
enough to choose Q 2 {z) = R(z), and the first-order phase integral approximation reduces to the WKB approximation. 
In the first-order approximation it is convenient to choose a root of Q 2 (z) as the lower integration limit in expression 
(fM)) . However, for higher orders, i.e. for 2N + 1 > 1, this is not possible because the function q(z) is singular at the 
zeros of Q 2 {z). In this case, it is convenient to express u>2 n (z) as a contour integral over a two-sheet Riemann surface 
where q(z) is single valued [17f . We define 



U2n(z) = 7; / Y 2n {z)Q(z)dz, 



(35) 



where t is a zero of Q 2 (z) and T t is an integration contour starting at the point corresponding to z over a Riemann 
sheet adjacent to the complex plane, and that encloses the point t, in the positive or negative sense and ends at the 
point z. 

If the function Q(z) is chosen conveniently, the quantity \i defined by 



\e(z)q(z)dz\ 



(36) 



is much smaller than 1. The function e(z) is given by the left side of Eq. (jT5J) 



eW _^.(„^ r V» W + ^_ 1 , (37, 

where the integral fi measures the accuracy of the phase- integral approximation (Tl| . 

We assume that the function Q 2 (z) is real over the real axis. Taking into account this restriction, we shall call 
turning point, the zero of Q 2 (z). We want to know the connection formulas at both sides of an isolated turning point 
z ret , i.e., a turning point which is located far from other turning points. We will adopt the terms "classically permitted 
region" and "classically forbidden region" in order to denote those ranges over the real axis where Q 2 (z) > and 
Q 2 {z) < 0, respectively. 

The connection formula for an approximate solution that crosses the turning point z ret from a classically permitted 
region to a classically forbidden region is [20( 
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q- 1/2 {z)\ cos (V(z)| + 1)^ \q-^ 2 (z)\ exp [\w(z)\ 



(38) 



The connection formula for an approximate solution that crosses the turning point z Te ± from a classically forbidden 
region to a classically permitted region is 



q~ l '\z)\ exp [- Kz)|] -> 2 |<T 1/2 (*)| cos (\to(z)\ - \) 



(39) 



It is important to emphasize the one-directional character of the connection formulas (|38p and (|39p . this means 
that the trace of the solution should be done in the direction indicated by the arrows in Eq. (|3"8"|) and Eq. {39} . 



A. The phase-integral power spectra 

We proceed to apply the phase integral approximation to Eq. ^ and Eq. Since the the square of the base 
function Q 2 S T (k, rf) possesses only one turning point rj s T that can be obtained after solving the equation Qg rj) = 
0, we can divide the axis rj into two regions: 



r) S T <-r]<0 Q ST < 0, 

-V < Vs.T, 



S.T > 0- 



The corresponding phase-integral solutions are: 
For rj ST < —rj < 



ut\k,ri) = ^ 



vt h% (k,rj) 



C-2 
dl 

2 

d-i 



q s 1/2 (k,rj) exp(- \w s (k,r))\) 

<?S 1/2 ( fc ' 7 ?) e xp(\ w s(k,r])\) i 
q T 1,2 {k,rj) exp (- \w T (k, T})\) 

— 1/2 

q T (k,r)) exp(\w T (k,r])\) . 



For -77 < r] S T 



u p k ht (k,rj) = ci 

+ c 2 

v P k h \k,ri) = di 

+ d 2 



<?S 1/2 ( fc ' 7 ?) C0S (\ W s(k, T))\ - j 

q s 1/2 (fc,?7) cos (\w s (k,T))\ + j 
q~ 1/2 (k,rj) cos (jw T (fc,?7)| - ^ 
q T 1/2 {k,r]) cos (\w T {k,r))\ + ^ 



(40) 
(41) 



(42) 



(43) 



(44) 



(45) 



where Ci, c 2 and di, di are obtained after comparing the asymptotic behavior of Eq. (|44[) and Eq. (|45[) with the 
asymptotic limit given by Eq. ([5]). In order to calculate the power spectra we substitute the growing part of the 
solutions (|42|) and (|43l into Eq. (J7]) and Eq. ((8]). We obtain the following expressions 



Pf\k) = Mm / '" ! N 



-fer^o 2tt 2 |z s (r/)| 5 



l^s*^ 7 ?)! exp (2 \ws{k,rj)\) 



pP hl (k) = lim 



fc 3 \<h\' 



k,^o 2tt 2 | (t7)|' 



|g T 1 (A:,77)|exp(2|w T (fc,?7)|). 



(46) 
(47) 



G 



The spectral indices ns and ut in the phase-integral approximation can be obtained respectively from Eq.([9]) and 
(fTOf and are given by 



nf\k) = 4+ lim 



n^ T hl {k) = 3 + lim 



d\n k 



din k 



din \q T 1 (k,r])\ 2 d\w T (k,i])\ 



dink 



dink 



(48) 
(49) 



III. APPLICATION TO POWER-LAW INFLATION 

The power-law inflationary model is a very simple model that allows one to solve the horizon and flatness problem. 
Since this model does not have a natural way of terminating the inflationary epoch, it is not physically acceptable, 
nevertheless its advantage lies in the possibility of analytically computing the solutions to the perturbation equations 
and the corresponding power spectra [2lL |22| . The power-law model allows testing approximations that are necessary 
in other inflationary models that do not exhibit analytic solutions. In this model, the scale factor is given by 



where v = | 
a > 0. 



a(t]) = lor)*'", (50) 
-^j . We have to impose the condition p > f in order that Eq. (f50"|) satisfies the inflationary condition 



Using the power-law scale factor (f50f we find that z$ = IqMpi 



Since, for this model, the differential 



equations governing the scalar and tensor perturbations are identical, we make the identification Uk = ffe with 



d 2 Uk 
drf 



k 2 - 



(v 2 - \) 



vr 



Uk = 0, 



(51) 



where the function Uk in Eq. (|5Tj) satisfies the boundary conditions ([5]) and ©■ 

Eq. (fBTj) can be exactly solved. The exact solution, satisfying the boundary conditions ([5]) and (JSJ) can be expressed 
in terms of a fractionary Hankel function [l5[ 



ufiv) = ^- ex P 
The exact power energy spectra are given by 

P§ X (k) 
P^ x {k) 

where 











-I) 











-nH<V{-k<n). 



„cx/ \ 1 3 — 2v 



7 2 M 2 f 
l IU Pl 



h ex (is)k 3 



(52) 

(53) 
(54) 



the corresponding spectral indices are 



3 e » 
h ex (y) 



l-2v 



r(iz) 



-1 2 



27rr(|) 



27rr(|) 



(55) 
(56) 



_2p_ 

p-V 



(57) 
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n§?(fc) = 2 - -i. (58) 
P- 1 

In order to apply the phase-integral method it is useful to introduce the variable z = krj. The function R(z) has the 
form 



R(z) = a + 

z z 



(59) 



where do = 1 and a_2 = | — v 2 are constants. In order to solve Eq. ([51]) with the help of the phase integral 
approximation we need to choose the base function Q(z). If we choose the square of the base function as Q 2 (z) — R{z) 
one obtains that the quantity fi(z, Zq) given by (|36[) is singular at the origin, which is the place where the boundary 
condition ^ has to be imposed. We can circumvent this problem making the following choice for the square of Q(z) 

El 



Q 2 (z) = & + 



(60) 



where bp and 6_2 are constants. The coefficients in Eq. (|60p are chosen in a way that makes the phase integral 
approximation valid to any order as — z — > 0. To verify the validity of the approximation at any order in a vicinity of 
zero, we require that the integral /i(z, zp) defined by Eq. (|36|) be finite as — z — > 0. For the problem we are discussing, 
it is enough to verify that the integral /i(z, zp) particularized to first order be valid in a vicinity of the origin. This 
condition guarantees that, at any order of approximation, the integral p,(z, zq) be defined in the vicinity of the origin. 
To first order, we have that q{z) = Q(z), and e = eo, then /i(z, Zp) can be written as 



H{z,zq) ■■ 

Making the corresponding computations we obtain 



\e (z)Q(z)dz\ 



(61) 



eo 



4a_ 2 - 46_2 - 1 



(46_ 2 - 4a_ 2 - 3) b a - b 



Ab 2 



z 2 + 0(z 4 ), 



(62) 



with 



4a_ 2 -4fe_ 2 -1 _i , 



Ah 



1/2 



(4fr_ 2 - 4a_ 2 - 7) bp a - bp 

b 1/2 

°-2 J 



z + 0(z 3 ). 



(63) 



With the help of Eq. (|6"3"|) and the expression (|6ip we get that the condition of validity of the first order approximation, 
as — z — > 0, implies that 



6-2 = a-2 - 



Using Eq. (|59 ]) -P? ]1 and Eq. (jM]) we have that 



(64) 



lim z 2 \Q 2 (z)-R{z)] =--. 
-2^0 J 4 



(65) 



The following election for Q 2 (z) is convergent as — z — > 0, and it is valid to any order of approximation 



Q 2 ( z )=R(z)-^. (66) 



The equation governing the modes k for the scalar and tensor is 

d 2 u k 



ij2 +Q z {z)u k = Q, (67) 
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where 



Q 2 (z) = l--z, (68) 

therefore, the phase-integral approximation is valid as — z — > 0, where the boundary condition (|6|) should be imposed. 

The square of the base function Q 2 (z) exhibits two turning points z re t = ±v. Since we are interested in the limit 
— z — > 0, we choose to work with the negative turning point. This turning point corresponds to the horizon k = aH 
(z = —v). The solution is defined in two ranges: On the left of the turning point, corresponding to scales lower than 
the horizon, we have the classically permitted region Q 2 (z) > where the solution oscillates. On the right of the 
turning point — v < —z < 0, corresponding to scales larger than the horizon, we have the classically forbidden region 
Q 2 (z) < 0, where the solution grows or decays exponentially. Fig. [T] (a) shows the two ranges where the solution is 
defined. In the phase-integral approximation, for —v < —z < 0, the solutions associated with the modes k for the 
scalar and tensor perturbations are given by and respectively. For — z < —v, the solutions are (gl]) and 





FIG. 1: (a) Behavior of the function Q 2 (z) for v < (b) Contour of integration Y- V (z) for — z < —v. (c) Contour of integration 
Y - v (z) for — v < —z < 0. The dashed line indicates the path of integration on the second Riemann sheet. 



Using Eq. ([32]) we have that the ninth-order approximation (2N + l = 9^7V = 4)of the function q(z) has the 
form 



4 

q(z) = J2 Y 2n Q{z) = (To + Y 2 + Y 4 + Y 6 + Y 8 ) Q(z). (69) 

In order to compute q(z) in Eq. (|6"9")l . we need to compute the coefficients Yz-n P6]) - ([3171) . The calculation of Y 2n 
requires the knowledge of the coefficients e„. Using Eq. (fl9|) and Eq. ([31) we obtain 



9 



z 2 (z 2 



z 1 (z 4 



4 V 2 

~2 



(z 2 — J/ 2 )' 



i- + 5^ a +2^ , (71) 



5|! + 28^ 2 z 4 + Mu 4 z 2 + 4i/») , (72) 

z 2 

- (6z 8 + 180t/ 2 z 6 + UOu 4 z 4 + 176v 6 z 2 + 8v 8 ) , (73) 

z 2 

g (30z 10 + 1320v 2 z 8 + 5400v 4 z 6 + 4576v 6 z 4 + 808v 8 z 2 + 16z/ 10 ) , (74) 

(z 2 — V 2 ) 
z 2 

w (I80z 12 + 10920i/V° + 67200j/ 4 2 8 + 98112i/ 6 z 6 + 38768z/ 8 z 4 + (75) 

(z 2 - j/ 2 ) 2 

3488j/ 10 z 2 + 32i/ 12 ) , 

z 2 

£ e = — ^TTI (l260z 14 + 100800^ 2 z 12 + 870240^ 4 z 10 + 1947456i/ 6 z 8 + (76) 



{z 2 -v 2 y ■ 

+ 1366416^ 8 z 6 + 291904^ 10 z 4 + 14560i/ 12 z 2 + 64z/ 14 ) . 
Inserting Eq ([70 |) -([7g ]) into Eq. J27|)-([HD we obtain that the first four 



z 2 



( z 2 _ v 2^ \\28 16 8 2 

z 2 /1073z 10 11171i/ 2 z 8 11511iAz 6 611i/ 6 z 4 213i/ 8 z 2 f 10 
( z 2 _ ^2)9 ^ 1024 + 256 + 64 + 4 + 8 + T~ 



_ z 2 /375733z 14 1752101z/V 2 7572705^ 4 z 10 4252813i/ 6 z 8 

8 ~~ ~ ( z 2 _ ^2)12 ^ 32768 + 2IM8 + 1024 + 256 + ^ ' 

1491943^ 8 z 6 39561^ 10 z 4 967^ 12 z 2 v 14 



128 16 8 2 

After computing the coefficients Y2 n up to N = 4 we obtain a ninth-order approximation for q{z). The next step 
is to compute to(z). Fig. [1] (b)-(c) show the contours of integration used for computing the integral u>(z) beginning 
with the second-order of approximation. Thus we have 

4 

Uj(z) = Wq{z) + 22 W 2 n(z), (81) 
n=l 

= f Q{z)dz + i V / Y 2n Q(z)dz. (82) 
J-u 2 n=i J r-» 

The expressions for u)2 n (z), written in the variable r), up to N — 4, are: 
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v arccos 



IW 



LU 2 

IUJ4 
iujQ 



= T 



(—krj) 2 
1 



v In 



—krj ' 



-krj < 



{-kr,f 



24 [±(-k V ) 2 T v 2 ?' 2 
1 

5760 [±{~-kri) 2 =p i/ 2 ] 9/2 
1 



[2i 



—kr\ 

3(-k v f] 



-v < —krj < 0, 



[375(-fcr ? ) 6 + 3654;/ 2 (-fc?7) 4 + 1512i/ 4 (-fcry) 2 - 16z^ 6 ] 



= T- 



322560 [±(-kr/) 2 =F f , 
-4744640i/ 4 (-fc?7) 6 + 1891200i/ 6 (-fcr ? ) 4 
1 



[67599(-fc?7) 10 + 1914210i/ 2 (-fcr ? ) 6 

21 / 



3440640 [±{-krf) 2 =F ^ 2 ] 21/2 
1965889800^ 4 (~fcry) 10 + 2884531440i/ 6 (-fcry) 8 



78720i/ 8 (-fc?7) 2 + 256/y 10 ] , 
[5635995(-/c?7) 14 + 318291750/v 2 (-fc?7) 12 - 
1135145088z/ 8 (-/b7) 6 - 



(83) 
(84) 
(85) 
(86) 
(87) 

(88) 



+99783936^ 10 (-/c?7) 4 + 881664t/ 12 (-/c7/) 2 - 2048i/ 4 ] , 

where the upper and lower expressions on the left side and the upper and lower signs on the right side in (I85|) - 
(|88l) correspond to —krj < —v and — v < —krj < respectively. After computing us(z), using the relations (|42|) 
and (|44[) we obtain a ninth-order phase integral approximation to the solution of the equation for scalar and tensor 
perturbations (15 1[) . The constants ci and ci are computed using the asymptotic behavior of Uk- We calculate the 

limit — krj — » 00 of the solution (f5Tjl on the left of the turning point (j44|) . Choosing C2 = — zci with ci = e l ( I/+5 ) ^ /V2 

and we obtain the asymptotic boundary condition given by Eq. ([5]). In order to compute the 

power spectrum we need to evaluate the limit — kr\ — > for the growing part of the solution (|42| . In this limit we 
have 



phi 



(v) 









exp 




3] 



where 



Using Eq. 



r-^(-* +i j;-i 



1 



1 



1260/v 5 1680i^ 7 

we have that the scalar and tensor spectra, given by Eq. (|4"6")) and Eq. (j4"Tj) are 

1 



where 



Pf\k) 
pP hl (k) 

9l hl = 
M hi = 



phi 7, 
,. ft, 



i 1.3-2;/ 



72 », f 2 6^ 
l O m Pl 

_UP hi h 3 ~ 2v 



l-2f 
3-2^ 

2 1 / 2 fP fl 

2^r~ 



jph'i 



-, 2 



2vr 



(89) 

(90) 

(91) 
(92) 

(93) 
(94) 



The spectral indices, given by Eq. (|4*5|) and Eq. (|4*9")) coincide with the exact ones given in Eq. ([57]) and Eq. ([55]) 
respectively. If we only keep the first term, — z/, in the exponential (|90[) . we obtain the first-order phase-integral 
approximation which coincides with the WKB method after using the Langer modification [1, 0, § ■ If we keep the two 
first terms in the exponential (|90p . we obtain the third-order phase integral approximation. It is worth mentioning 
that, for the power-law model, the tensor and scalar spectral indices do not depend on the order of approximation. 
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IV. RESULTS 



In this section we proceed to compare the analytic solutions for the wave function Uk and the scalar and power 
spectra with the results obtained using the phase integral approximation. Fig. [2] compares the real part of the 
analytic solution of Uk with the ninth-order phase integral approximation. The plot is made using the number of 
e-folds, N = In a, as the independent variable. As expected, the phase integral approximation solution diverges at 
the root of q(rj) 



0.00 



CD 

DC 



-2.00 



-4.00 



-6.00 






e-folds 



FIG. 2: ?R(uk) for the power-law inflationary model with p — 10 with k — 1.389ftMpc . Solid line: analytic solution; dot-dashed 
line: ninth-order phase integral approximation. 

Now, we compare the scalar and tensor power spectra calculated using the phase-integral approximation with the 
results obtained with the slow-roll and uniform approximation methods. From Ref. t 8j (Eq. (63) and Eq. (64)) we 
obtain that the scalar and tensor power spectra in the slow-roll approximation are 

p l r ( k ) = l2^FT9t r k 3 ~ 2 », (95) 



llM* Pl 
1 



(96) 



with 



g s J = [1 + 2(2 - In 2 + 6)(2e + S) - 2e] 



l-2v 
3-2u 



2~ v \l-2v 



1/-1/2 



2tt 



h s J = [l-2(ln2 + 6-l)e] 



2 1 /2-"|i_2i/| 
2^ 



I/-1/2 



(97) 



(98) 



where b is the Euler constant, 2 - In 2 - b ~ 0.7296, In 2 + b - 1 ~ 0.2704, and e = -S = i. In the slow-roll 
approximation e ^ 1, therefore, for the power-law model, the slow-roll approximation is better suited for large values 
of the parameter p. 
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Using the result obtained in Ref. [6|] (Eq. (109)), we obtain an expression for the second-order uniform approxima- 
tion for the power spectrum associated with the scalar and tensor perturbations. They have the form: 



PT(k) 



I 2 M 2 

l O m Pl 



ua L,3-2u 
9v ft ) 



ua 7,3 — 2^ 



(99) 
(100) 



with 



1+ 6u) \ 3-2v 



hT = 1 



(21/) 



I/-1/2 p -v 



2tt 



(101) 



2tt 



Omitting the factor 1 /Qv in Eq. IjlOip we obtain the first-order uniform approximation, result that coincides with the 
first-order phase- integral approximation and the WKB method with the Langer modification Q . Keeping the second 
term of the expressions in Eq. (|101|) one gets the second-order uniform approximation. 

We want to compare the analytic expression for the scalar and tensor power spectra for different values of k with 
the ninth-order phase integral approximation, the slow-roll approximation and the first and second-order uniform 
approximation. Fig. [3] and Fig. [4] show the power spectra Ps(k) and Pr(k) calculated analytically as well as with 

different approximation methods. The inset in each figure represents an enlargement that permits one to compare 

tl .. . .. . - „, 

P 




FIG. 3: (a) Scalar power spectrum Ps(k) and (b) tensor power spectrum Pr(k) for the power-law inflationary model with 
p = 10. The solid line indicates the analytic solution. The dot-dashed line: ninth-order phase integral approximation. The 
inset in each figure represents an enlargement that permits one to compare the accuracy of each method. 



Fig. [5] shows the quotient between the exact analytic result and the result obtained using different methods of 
approximation g 3 p /gp X for different values of p. Using the ninth-order phase integral approximation we obtain a 
horizontal line equal to unity. We observe analogous behavior when we consider the quotient h) p /h c *. As it was 
already mentioned, the slow-roll approximation is better suited for large values of p because e = 1/p. The WKB 
method gives a better approximation for small values of p because the condition /i <C 1 is valid in this case. It should 
be noticed that the first-order phase-integral approximation, the WKB method with the Langer modification and the 
first-order uniform approximation give the same result. 
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FIG. 4: (a) Scalar power spectrum Ps (k) and (b) tensor power spectrum Pr(k) for the power-law inflationary model with p — 10. 
Solid line: Analytic solution; dotted line: slow-roll approximation; dashed line: first-order phase-integral approximation, WKB 
and first-order uniform approximation; dot-dashed line: third-order phase integral approximation; dot-dot-dashed line: second- 
order uniform approximation. The inset in each figure represents an enlargement that permits one to compare the accuracy of 
each method. 



Fig 6(a) and Fig |6(b)| show, the relative error of Ps(k) and Pr(k) respectively for different approximation methods. 
The relative error is computed using the following expression: 



rel. error Ps,t(^) 



papprox 



(k) 



Psj'(k) 



P^ T (k) 



100. 



(102) 



Approximation 


|rel. error Ps(A;)| 


|rel. error Pr(k)\ 


phi9 a 


0.0015% 


0.0015% 


phi3 b 


0.1% 


0.1% 


ua2° 


0.4% 


0.4% 


sr d 


0.5% 


0.5% 


phil°, WKBU, ual g 


10% 


10% 



TABLE I: Relative error obtained using different methods of approximation for the power-law inflationary model with p — 10 
for the mode k = 1.369 ftMpc" 1 . 

"Ninth-order phase integral approximation. 

b Third-order phase integral approximation. 

°Second-order uniform approximation. 

d Slow-roll approximation. 

Tirst-order phase integral approximation. 

f WKB approximation with the Langer modification. 

Tirst-order uniform approximation. 



Table [J shows that the ninth-order phase integral approximation gives the smallest relative error for Ps(k) and 
.Pr(fc). The third-order phase integral approximation gives a relative error smaller than the one obtained using the 
second-order uniform approximation. 

Finally, some words about the accuracy of the phase-integral approximation are in order: The smallness of the 
integral fj,(z, Zq) given by Eq. (|36|) is a measure of the accuracy of the phase-integral approximation. Fig. [5] shows 
that the phase-integral method fails in the vicinity of the turning point —v, range where the /^-integral diverges. The 
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10 11 12 



FIG. 5: Evolution of the ratio g J (p)/g ex (p) versus p. Solid line: ninth-order phase integral approximation; dotted line: slow- 
roll approximation; dashed line: first-order phase integral, WKB and first-order uniform approximations; dot-dashed line: 
third-order phase-integral approximation; dot-dot-dashed line: second-order uniform approximation. 



selection of the base function Q(z) given by Eq. ([55]) guarantees that far from the turning point at any order 

of approximation. Since the scalar and tensor power spectra as well as the spectral indices are evaluated as — kr/^ 0, 
the limit is taken far from the horizon (turning point), therefore their computation is not affected by the presence of 
the turning point. For p = 10 and k = 1.369/iMpc - 1 the turning point is v = — 1.61Mp, , and the power spectrum is 
evaluated at krj = — 10 _8 Mp i 1 , a point in the limit where the function Uk exhibits the asymptotic behavior ((6|). Fig|3j 
Fig. Q] and Table U show the accuracy of the phase- integral approximation. 

In the present paper we have shown that, in comparison with other approximation methods, the phase integral 
approach gives very good results for the scalar and tensor spectra in the power-law inflationary model. We have also 
seen that the phase integral approach reproduces the exact spectral indices in the power-law model. Since the WKB 
method can be regarded as a first-order approximation of the phase-integral approximation with Q 2 (z) = R(z), it 
should be expected that the phase-integral method works in those cases where the WKB methods gives good estimates 
and slow-roll fails, that is the case where inflation is generated by a chaotic potential with a step [23|]. The good 
agreement between the analytic results and the phase-integral approximation shows that the phase integral method 
is a very useful tool for computing the scalar and tensor power spectra in more realistic inflationary scenarios, in a 
forthcoming work we will implement the phase integral method to quadratic and quartic chaotic models. 
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